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Abstract 

Among the problems still open in the study of stellar 
structure, we discuss in particular some issues related to 
the study of convection. We have recently built up com- 
plete stellar models, adopting a consistent formulation of 
convection both in the non gray atmosphere and in the 
interior, to be used for non adiabatic pulsational analy- 
sis, and discuss in some details two problems which have 
been clarified by these models: the physical interpreta- 
tion of one of the main sequence "Bohm-Vitense gaps", 
and the necessity of parametrizing pre main sequence con- 
vection differently from the main sequence convection. 
We also report preliminary results of the application to 
the solar model of the non-local turbulence equations by 
|Canuto fc Dubovikov (1998) in the down-gradient approx- 
imation. 
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1. Introduction 

The list of unsolved problems in stellar structure is long: 
microscopic diffusion; mass loss; rotation, its evolution 
and its effects on chemical mixing; magnetic fields, and 
their interaction with rotation and convection; and an ad- 
equate model of convection able to describe the overa- 
diabaticity in the shallow layers, and non-local aspects 
like overshooting. Asteroseismology, and especially aster- 
oseismology from the space, is destined to be very useful 
in providing us with additional constraints to our mod- 
els. Here, we only touch the problem of convection in 
stars, and in particular the problem of envelope convec- 
tion. The convective envelope in low mass stars has a fun- 
damental role in <5Scuti, 7Doradus and solar-type oscil- 
lations. Furthermore, the boundaries between the corre- 
sponding domains of instability in the HR diagram seem 
to be linked to changes in the convection features. The 
non adiabatic seismological analysis requires knowledge 
of the entire stellar structure up to the atmosphere, and 
selfconsistent convection models for the atmosphere and 
interior integration are necessary for this use. On the other 
hand, progress in the study of convection is slow enough 
that today's main attitude with respect to this problem 



is the following: parametrize in a rough way the overa- 
diabatic convection, by using constraints either coming 
from the analysis of stellar spectra or from other struc- 
tural observational information (e.g. the p-modes solar 
patterns, or the depletion of light elements in the stel- 
lar envelopes, or empirical radii constraints). This can 
be sufficient for studying solar oscillations, but it is dif- 
ficult to extend models so much parametric to regions of 
the HR diagram in which the observational constraints 
are scarce (or none). We have been attempting to build 
up complete stellar models to be used for non adiabatic 
pulsational analysis. These models have clarified some in- 
teresting problems in stellar evolution which we discuss 
here in more detail: i): they have provided an interpreta- 
tion for the "Bohm Vitense gap" observed in open clus- 
ters at B-V-0.35 IjRachford fe Canterna 2000|) : it): they 
have made explicitly clear the necessity of parametrizing 
pre main sequence convection differently from the main 
sequence convection, possibly implying that there is an 
hidden "second parameter" acting to change the convec- 
tion behaviour during the first phases of stellar evolution. 
Finally, we show the results of preliminary computation 
of non local convection in the solar model, which indicate 
that the overshooting expected at the bottom of the solar 
convective envelope is ~ 0.02Hp. 

2. Convection in complex evolutionary phases 

A good example of complex models of stellar structure are 
the intermediate mass Asymptotic Giant Branch (AGB) 
stars, those which evolve through the phase called Hot 
Bottom Burning (HBB). These stars are extremely impor- 
tant for the chemical evolution of galaxies, and they are 
probably the key ingredient to develop primordial chemi- 
cal inhomogeneities among Globular Clusters stars. Dur- 
ing 90% of their life, they are fueled by hydrogen, whose 
stationary burning in an external shell is interrupted by 
the sudden ignition of the helium beneath. During the en- 
suing thermal runaway (the Thermal Pulse -TP- phase), 
the hydrogen envelope expands, hydrogen burning stops 
and external convection reaches into the helium nuclearly 
processed layer, provoking the 'third dredge up' phase by 
which nuclearly processed material appears at the stellar 
surface. In the most massive AGB stars (4-8M Q ), the 
stationary hydrogen burning shell is partially contained 
into the convective envelope, that is, the bottom of the 
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convective envelope attains temperatures so large that nu- 
clear reactions take place there. The CN -or even CNO, 
for low metallicity stars IjVentura et al. 2001J) - nuclearly 
processed material is convected to the surface and then 
convected back to the bottom. So the nuclear processing 
which can be observed in the atmospheres of these stars is 
due to the combined action of HBB and of the third dredge 
up. In spite of the fundamental importance of these ob- 
jects evolution, for sure this phase is not easily modeled! 
One of the big problems is the (unknown) mass loss, which 
affects the yield of the processed chemistry both through 
the time dependent loss of the envelope into the interstel- 
lar medium, and through altering the temperature struc- 
ture of the envelope. In addition, these stellar models are 
heavily dependent on the way we model convection it all 
its aspects: 

1. the convective model, that is the convective flux and 
temperature gradients computation; 

2. overshooting, that is mixing outside the formal con- 
vection boundaries. For AGBs, this includes both the 
problem of core-overshooting, which affects the rela- 
tion between initial mass and mass at the beginning of 
the TP phase, and the possible overshooting below the 
convective envelope, which affects the nucleosynthesis; 

3. time dependence: mixing must be treated as non in- 
stantaneous, and coupled with nucleosynthesis, for all 
the elements for which the convective mixing timescale 
is of the same order of the nuclear burning timescale. 

3. Convection in more "normal" stars 

In the present discussion, we leave entirely aside the "big 
problems" of convection modelling — apart from the last 
section. While turbulence in stars is compressible and non 
local, we must accept that for a long time, for general pur- 
poses of computing stellar structure for any mass, chem- 
istry and evolutionary phase, we still will deal with incom- 
pressible and local models. 3D radiation hydrodynamics 
(RHD) simulations still miss the computer power needed 
to deal with deep envelope convection, although great in- 
sight has been obtained in the atmospheric studies. In- 
teresting results are available by 2D simulations. Analytic 
non-local convection models have recently been applied to 
stellar atmospheres of A stars ( |Kupka fc Montgomery 20021 
and we may foresee important developments also in this 
area. Among the local models, the Mixing Length The- 




ory (MLT) by Bohm-Vitense (1958) certainly is still dom- 
inant, in spite of its limitations. Other local models have 
become available, among them the Full Spectrum Turbu- 
lence (FS T) model by|Canuto fc Mazzitelli (1991)| and its 
variant in Canuto et al. (1996)[ which have been adopted 
for a variety of applications in the latest 10 years. 

As no available theory is self consistent, the domi- 
nant attitude in stellar studies -unless there are particular 
requirements- is to use the convective model as a 'black 
box' simply to infer the stellar properties in the region in 



Figure 1. Stratification of temperature vs. pressure in the 
atmosphere and sub-atmosphere of solar models by Mon- 
talbdn et al. (2003). Track 1 is the model employing FST 
convection both in the atmosphere and interior, the match 
is done at T pn = 10; models from 2 to 4 have MLT at- 
mospheres with a atm = 0.5, and MLT interiors with «i n 
chosen so to fit the solar radius. The differences are due 
to the different choices of matching points: 2: r pn = 100 
and ain = 6.3; 3: tau pn = 10 and a- m = 2.3; 4 : r ph = 1 
and a; n = 1-75. Model 5 employes the AH97 atmospheres 
down to T p h = 100, and has a- m = 1.9 (the model shown 
in Figure^). We see that the subatmo spheric structure is 
very different for the different models, in spite of the 'solar 
radius' calibration. 



which convection becomes adiabatic. In general, in fact, 
the main important property which a convection model 
is required to provide is the stellar radius: whatever is 
in the 'black box' of convection, we first of all wish to 
know the radius (or the T e g, or, equivalently, the entropy 
jump between the adiabatic interior and the surface). The 
solar radius in particular is generally adjusted in the mod- 
els by varying the convective efficiency. In the MLT this 
can be done by varying the 'mixing length' parameter, 
that is the ratio of mixing length to pressure scale height 
a = I /H p . Recently it has been stressed the importance 
of surpassing the Eddington approximation, or the gray 
atmosphere approximations to correctly describe the op- 
tical atmosphere. For many stellar situations, a non gray 
model atmosphere — which also must reproduce the ob- 
served spectral features — includes as necessary ingredient 
a treatment of convection, which, in the MLT, will be char- 
acterized by a given a a t m down to the optical depth r p h at 
which the atmosphere will be matched to the interior com- 
putation. In the interior, convection will be characterized 
by a given a- m . Thus, actually, a 'modern' MLT structure 
will have three parameters: a a tm, cc ln and T pn . Figure 
illustrates this classic problem by using our recent models 
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(Moutal ban et al. 2003(1 which adopt as boundary condi- 
tions the NextGen atmospheric grid IjAllard et al. 19971 
hereinafter referred to as AH97), computed by assuming 
MLT convection with a a t m = 1, down to optical depth 
Tph = 100. 

Adopting ct; n = 1.9 for the interior computations, we 
obtain a set of tracks, among which the solar track passes 
through the solar location at the solar age. The solar mass 
track with a- m — 1.0, on the contrary, is ~ 380K cooler 
at the solar luminosity. This well known fact remembers 
how large is the variation in T e g allowed by changing the 
MLT parameter. The meaning of this solar radius adjust- 
ment is simply the following: when using the MLT, the 
entropy jump, necessary to fit the solar T e ff, between the 
adiabatic layers and the surface corresponds to what is 
obtained with the MLT and a = 1.9 up to r p h = 100, 
and with MLT and a — 1 up to the surface. Of course, 
the model can not tell us anything about the temperature 
gradient layer by layer within the envelope. If we only 
wish to use convection as a 'black box', we just take the 
a- ln — 1.9 model for the Sun 1 . Of course, our assumption 
will only be valid for the Sun itself. 

An interesting broadening of of this point of view to a 



technically necessary for stellar stability studies, to avoid 
discontinuities in the physical quantities. A prototype of 



these models has been discussed by Schlattl et al. (1997) 



wider part of the HR diagram has been given by Ludwig et 



who performed detailed 2D numerical RHD calculations 
of time-dependent compressible convection, in the range 
4300< T cff < 7100K, 2.54< \ogg <4.74 (cgs) for solar 
composition. They used these models to 'calibrate' the ef- 
fective a for these envelopes, that is the value of a which 
provides the same specific entropy jump (using a gray at- 
mosphere, therefore this approximation defines a unique 
average a parameter) than their 2D models. They found 
values from 1.3H p for F dwarfs up to 1.8H p for K giants. 
This calibration of a, again, does not tell us anything 
about the temperature gradient layer by layer. 

The approach of keeping convection as an entirely black 
box is very useful, but it has -obviously- many limita- 
tions. In particular, if we wish to study the excitation 
of oscillation instabilities in stars we must know the full 
stratification of physical quantities in the star, up to the 
atmosphere. 

4. Convection modelling in full stellar models 

We have to recognize that no general purpose convection 
model is presently available, which can be meaningfully 
applied to any stellar structure in order to know the layer 
by layer stratification of physical quantities. Nevertheless, 
recently there have been several attempts to produce en- 
tire stellar models, which satisfy some observational con- 
straints on the atmospheric and envelope convection, and 
in which the atmospheric and the envelope convection are 
matched 'smoothly' each other. A smooth match is at least 

1 and generally we forget, or at least do not even mention, 
that a is actually a a tm = 1 for most of the overadiabatic region, 
which is contained in the model atmosphere grid. 



for the Sun. They built up a non gray ID model atmo- 
sphere with aatm = 0.5 down to r p h = 20, and matched 
it to a subatmosphere and MLT interior, guided by the 
results of 2D hydro models. To do this, they had to vary 
the ain parameter in the interior computation, with the 
aim to provide the solar radius and obtain a temperature 
stratification similar to that of the 2D models. The aim 
of this study was to explore the influence of the physical 
inputs on the solar p-modes. 

Of course, for the Sun we have an enormous number of 
constraints — from the p-modes themselves, and from the 
precise knowledge of the solar radius — which help to build 
up a fully parametric model. But what should we do to 
extend the analysis to other stars? What kind of model at- 
mosphere can we use, and how do we produce a 'meaning- 
ful' match of the atmospheric and interior computation? 
There are not yet many model atmosphere grids available, 
and many of them have been computed for "atmospheric 
purposes" only, that is to produce an adequate modelling 
of the re gion from which the stellar spectum emerges, so 
a l-tl(i€^99} c not entirely apt to be used as boundary con- 
ditions for full stellar models. For instance, the quoted 
AH97 models, computed with the PHOENIX code, adopt 
MLT convection with a = 1 and a total of 50 layers down 
to t == 100. Heiter et al. (2002) have adopted a new ver- 
sion of Kurucz's <|Kurucz 1993|l model atmospheres (the 



Vienna- ATLAS9 code), in which the convection model is 
either the MLT with a = 0.5 or the FST (both in the 
ICanuto &: Mazzitelli 199T1 and in the ICanuto et al. 19961 
versions), and they have increased the number of layers 
in the latter models to 288. We have recently used all 
these grids of model atmospheres to explore the meaning 
of the match between atmosphere and interior, and what 
happens when one extends to other models (pre main se- 
quence, and giants) the approximations adopted for the 
solar model. The extensive results of this study are pre- 
sented elsewhere (Montal ban et al. 2003|l . Here we show, 
in Figure ^ that the different solar models which can be 
built up, all satisfying the constraint of the solar fit, can 
have very different subatmospheric structures, depending 
on the choice of the convection efficiency in the different 
layers. 

In Figure \l\we also see that the model adopting the 
FST convection both in the atmosphere and in the inte- 
rior does not show discontinuities in the T vs. P stratifi- 
cation. The structure smoothly passes from a very ineffi- 
cient convection in the atmosphere, similar to that of the 
MLT models with a = 0.5, to such an efficient convec- 
tion inside, that the precise solar fit is easily obtained (see 
e.g. ICanuto fc Mazzitelli 199T)l 2 . Computing "full FST" 



2 Remember that, whatever the choice of the free param- 
eters in the FST convection, it is not possible to obtain a 



4 




Figure 2. Temperature vs. pressure stratification in mod- 
els from 2 to -ZMq at 10 8 yr, built with the Vienna- 
ATLAS9 FST boundary conditions. Both the atmosphere 
(dot-dashed) and the interior (full lines) are shown. The 
transition between almost radiative envelope models ( down 
to J.^Mq ) and models with well developed envelope con- 
vection (below 1.42Mq ) is clearly shown. The same tran- 
sition is much smoother in models with MLT convection, 
having also in atmosphere a larger than 1. 



Figure 3. HR diagram of the Hyades from Hipparcos data 
(Debrujine et al. 2000, 2001), and an isochrone of 600Myr 
built by the present models (Montalbdn et al. 2003) com- 
pared with an isochrone built by gray models (Ventura 
et al. 1998a). MLT models with a = 1.6 which fits the 
Solar location have a smoother transition (Mazzitelli & 
D'Antona 1993). We also show the 'lithium dip' (Boes- 
gaard & Tripicco 1986) (triangles), and the logarithm of 
the convective mass (scale to the right). 



structures, then, has as a first a "formal" advantage, as it 
helps to avoid the problem of discontinuities in the physi- 
cal quantities for instability studies. In addition, there are 
other, more physical, reasons to use this model as probe: 
a series of works 011 1) T e R determination from H Q and 
lines fe.g. lFiihrmann et al. 1993llvan't Veer-Menneret et al. 
2): theoretical predictions of H Q and from ID models 



the HR diagram. We give here one positive example (the 
interpretation of the Bohm-Vitense gap at B-V~ 0.35) 
and one "negative" example, namely the impossibility of 
explaining the patterns of Lithium depletion in the Sun 
and in open clusters with the FST model. This latter re- 
however, mav be telling us something else on the 
behaviour of convection in young convective stars. 



i|Gardiner et al. 1999jl and from 2D models ( |Steffen &: Ludwig 1 999); 



3): theoretical predictions oib — y and c Stromgren indices 
(Smalley & Kupka 1997) and abundance determinations 
l|Heiter et al. 19 98) indicate that, even if a ID, homoge- 
neous model can not explain all the spectroscopic and pho- 
tometric observations, model atmospheres which predict 
temperature gradients closer to the radiative gradient are 
in better overall agreement with the observations. These 
arguments favors those models in which convection in the 
atmosphere is less efficient than predicted by models hav- 
ing a > 1), and led Heiter et al. (2002) to adopt either 
the FST model or the a = 0.5 MLT models to compute 
their atmospheric grids. These were also further motiva- 
tions to produce "full FST" stellar models by use of their 
atmospheric grids. 

We may ask whether the "full FST" models have shown 
features interesting enough to render it useful to explore 
the application of the FST convection in different parts of 



lM© location, at L©, with T e ft much smaller or larger than 
the solar T e g, contrary to the MLT models. 



5. The Bohm Vitense gap at B-V~ 0.35 

One of the most interesting results in our recent explo- 
ration of "full" FST models regards the main sequence: the 
FST convection, due to its very low efficiency close to the 
stellar surface, and high efficiency in the inner subatmo- 
spheres, yields a very sharp transition between structures 
which are convective only in the surface layers, and struc- 
tures which show a well developped convection also in the 
interior (D'A ntona et al. 2 002) . This is shown in Figurc[2] 
by comparing the stratification of models of different mass 
in the P-T plane: models down to 1.42Mq are mostly ra- 
diative (or convection is so inefficient that the convective 
gradient sticks very close to the radiative one), while sud- 
denly, at 1.41M , an extended adiabatic convection re- 
gion appears for a larger part of the envelope. The fast in- 
crease in the convective mass fraction (Af ce ) as a function 
of the main sequence color B-V is shown in Figure 3. This 
characteristic of the models is reflected in their mass-T e g 
relation, which becomes suddenly steeper around T e fj~ 
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Figure 4- The Hyades Hipparcos HR diagram is compared 
with two simulations containing 100 stars between 1.2 and 
1.7Mq , randomly extracted from a Salpeter mass func- 
tion distribution, following the FST or MLT mass-color 
isochrone of 600 Myr. The MLT simulation is shifted by 
+0.7 mag, and the FST one by 1.2 mag for clarity. At 
the bottom, the FST (dashed) and MLT color histograms 
are shown (scale at the right). The FST mass-color dis- 
tribution produces a 'gap' at B-V^0.35, while the MLT 
isochrone does not. The good correspondence with the 
Hyades gap indicates that the transition between structures 
which are convective only in the surface layers, and struc- 
tures which show a well developped convection also in the 
interior is very sharp 



Figure 5. The general HR diagram of the computed tracks, 
which extend from the PMS to the MS and then to the 
Giant Branch. The tracks of 1.2, 1.4, 1-45, 1.5, 1.55, 1.6, 
1.65M.Q are shown. The diagonal (green) line on the right 
shows schematically the separation line between higly con- 
vective structures on the right and structures with convec- 
tion limited to the atmosphere (on the left). The diago- 
nal dotted (black) line indicates the observational red edge 
of the 5 Scuti instability strip, according to Pamyatnykh 
(2000). The two dash-dotted (red) lines indicate the obser- 
vational boundaries of the 7 Dor strip (Zerbi 2001). The 
MS at lCPyr age is also shown. 



be very different for stars having so different convective 
structure. In particular, it could represent the dividing line 
between coherent pulsations and solar-type oscillations. In 
this case, we should have expected that also the stars be- 
tween the red edge of the 7 Doradus instability strip and 
the transition line are pulsating. Notice that the location 
of the red edge of the 7 Doradus strip is based on observa- 
tions from the ground and might still be uncertain in the 



6800K, and is also apparent in the HR diagram as a sud- 
den change of slope (Figure 3). Using numerical simula- 
tions we have shown that this feature produces a stel- 
lar depletion which is consistent with the gap seen in the 
Hyades at 0.33 £ B-V ^ 0.38, one of the so called "Bohm 
Vitense gaps" after Bohm-Vitense & Canterna (1974) and 
Bohm-Vitense (1982)| found by |Rachford fc Canterna (2000) heoretical HR diagram. On the other hand, many other 

structural parameters may have a role in defining these 
instability strips, and at least two of these — elements dif- 
fusion and rotation — should be considered. A further note 
of caution on the naivety of this proposal comes also from 
the complex behavior of the chromospheric and transition 
layer indicators for the MS stars on the right of the tran- 
sition line IjBohm- Vitense et al. 2002 ) . 

Will the Bohm Vitense gap appear also in the field 
stars? Here the problem will be much more complicated by 
the blurring due to the fact that we have a mixture of ages 
and metallicities! Figure[7]shows simulations for three stel- 
lar populations with different metallicities (Z=0.01, 0.02 
and 0.03), covering ages from 10 7 to 10 9 yr, and distributed 
according to a Salpeter mass function. The diagonal lines 
separate the location of the models having deep envelope 



in 6 our of 9 open clusters which have been investigated. 
The standard MLT models do not show this behavior (see 
Figure HJ. 

The very sharp variation of the stellar structures in 
the HR diagram (or in the T G ff gravity plane) is shown as 
a transition line in the HR diagram of Figure Models 
on the right of this line have extended convective regions, 
and models on the left have very small convective regions, 
independently from their evolutionary phase (pre, on or 
post the main sequence). This can be seen in Figure^] The 
transition line is compared in Figure [5] with the location 
of the S Scuti and 7 Doradus instability strips. We may 
speculate that our transition line separates HR diagram 
regions which harbor different modalities of stellar oscilla- 
tion patterns, as the excitation or driving mechanisms can 
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Figure 6. The radiative (black), adiabatic (red) and supera- 
diabatic (blue) gradient V — V a d in the external layers of 
the selected models along the 1.5M@ evolution (shown 
as open squares in Figure [3| versus the logarithm of the 
pressure. The first and third model have their convection 
region limited to the atmospheric layers, while the second 
one shows an extension of convection down to log P ~ 9, 
well inside the envelope. The transition occurs for a very 
small variation of the physical parameters, and is mainly 
due to the high inefficiency of the FST model in the at- 
mospheric regions. 



However, a different problem emerges from the most 
recent computation of solar models: they deplete too much 
lithium during the PMS evolution l|D'Antona fc Mazzitelli 19941 
IIVAntona fc Mazzitelli 1997lEfrhlat,t,l fa WeisH999llPia,u fa Tnrcl 
and are incompatible with the young open clusters ob- 
servations. This problem is most severe in models using 
very efficient convection models, in fact it is more relevant 
for the D'Antona and Mazzitelli (1994 and 1997) models 
adopting the FST convection. MLT models of the most re- 
cent generation, adopting updated equations of state and 
opacities also deplete too much lithium. D'Antona & Montalban (z 
have recently reappreciated that the problem is severe in 
tracks whose convection is adjusted to provide the solar ra- 
dius at the solar age! If one does not require the solar fit, it 
is easy to decrease the convection efficiency and to obtain 
a negligible PMS lithium depletion. In Fig.[S]we compare 
the lithium depletion predicted by many of our models, 
computed using AH97 atmospheres (a atm = 1) as bound- 



ary conditions, with the Pleiades data by |Soderblom et al. (1993) 
and 



Garcia Lopez et al. (1994) Only the upper squares, 



corresponding to the models employing a ln = 1 , which do 
not fit the Sun, are compatible with the data. The HR 
diagram of the models is shown in Figure [5] The a = 1 
tracks are on the right, much cooler than what is needed to 
allow the solar fit! In addition, Figure|H|opens up the com- 
plex problem of the location in the HR diagram of PMS 
tracks (for a review see, e.g. ID'Antona et al. 20001 while 
our full recent computation for this phase can be found in 
Montalb an et al. 2003|l . When we compare the location of 
PMS theoretical tracks with the observed few data of PMS 
stars for which an independent determination of mass is 
available 3 , the tracks most consistent with the observa- 
tions are again those with cooler atmospheres (higher mass 
for a given spectral type) and thus those which, generally, 
provide a radius larger than R Q for the solar model. In 
fact, Figure shows that the a — 1 tracks are well com- 
patible with the location of the secondary component of 
the eclipsing spectroscopic binary RXJ 0529.4+0041, one 
of the best determined PMS masses l|Covino et al. 20011 
It is then clear that 1 ) the HR diagram location of the 
tracks during the PMS evolution, and 2) PMS lithium de- 
pletion, are two problems correlated each other, as it could 
have been expected, because, the smaller the T e s of the 
Hayashi track, the smaller is the temperature at the base 
of the convective envelope during the possible phase of 
lithium burning. Both the HR diagram location and the 
system abundance (log N(Li)=3.31±0.04 l|Anders fa Gr cvcsseiMTO ji. depletion seem to be compatible only with mod- 



convection (at cooler T e g) from those, hotter, having very 
thin atmospheric convection. A gap is apparent on the 
left of each line, which shift to cooler T e g with increasing 
metallicity, so that it will not be easy to locate a gap in a 
sample of stars not homogeneous in metallicity. 



6. The 'historical' problem of Lithium 

In the years 1965-1990, we called 'problem of lithium in 
the Sun' the inhability of the solar mass evolutionary 
tracks to burn a substantial fraction of their initial lithium 
during the Pre Main Sequence (PMS). This result was gen- 
erally taken as a good proof that additional mechanisms 
for depletion were required, acting during the long solar 
MS lifetime, to reduce by a factor ~140 the initial solar 



This interpretation still today is taken as most plausible 
one, confirmed by the scarce lithium depletion, at the solar 
mass, in young open clusters (see e.g. |Chaboyer 1998| ). In 
fact the lithium vs. T c s relation for the MS stars of young 
open clusters indicates a lithium depletion by at most a 
factor two for the solar mass in young clusters, while it is 
compatible with the solar depletion in some stars of the 
cluster M67, close to the solar age. For recent reviews see 
e.g. 



els in which PMS convection is much less efficient than 
MS convection. Is this simply another proof that we are 
not able to model convection, or that there are unsolved 
problems with the opacities? Or does it mean that there is 



Jeffries (2000) and Pasquini (2000) 



3 these stars either belong to binaries (e.g. 
ICovino et al. 200T1 |Steffen et al. 2001)| or one can mea- 
sure the stellar mass from the dynamical properties of their 
protoplanetary disks (Simon et al. 200011 . 
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some other parameter playing a role in the PMS -and not 
on the MS? It is probably too early to derive strong con- 
clusions, but we have suggested that PMS convection is 
inhibited by the thermal role of a dynamo built magnetic 
field l| Ventura et al. 1998bl ID'Antona et al. 2000jl . 

7. Non-local convection in the Sun 

Always in the framework of stellar evolution for astero- 
seismology, we performed one more consistency test, of a 
completely different nature. In this case, our attention is 
not at the surface, but at the bottom of the solar con- 
vective zone (CZ). In fact, a constraint to the thickness 
of the overshooting layer in this region has been set by 
helioseismology (e.g. IBasu fc Antia 1997|1 . and it can not 
be larger than ^0.05_ff p . We constructed a detailed FST 
solar model matching the correct thickness of the CZ (the 
surface boundary condictions are in this case of negligi- 
ble importance) and applied the treatment suggested by 
ICanuto fc Dubovikov 1 998 (CD98) to a thin region cen- 
tered around the formal Schwarzschild boundary of con- 
vection, to gain insight on what happens when a fully 
non-local turbulence theory is applied to a stellar struc- 
ture. 

More in detail, we have computed from the local model 
the starting distribution for the quantities: K (turbulent 

2 

kinetic energy in the radial direction), 9 (mean quadra- 
tic temperature variance), J (convective flux), w (mean 
quadratic velocity variance) and e (dissipation rate), ac- 
cording to (42a-c), (43a-c) and (44a-d) in CD98. Then, 
temporal relaxation to the five above quantities has been 
allowed, according to the equations (19a-d) and (35a-b), 
until stationary conditions were reached. For each relax- 
ation step, the gradient /?: 



^ = ^(V-V od ) 
tip 

has been updated according to: 
J + F(K) 



P = Prad 



X 



(1) 



(2) 



where x is the thermometric conductivity, ft ra d — T(V ra( 2— 
V 'ad) I 'H p . Further, F(K) is the kinetic energy flow ~^^r 

K 2 F 

and v t is the turbulent viscosity, v t oc — . 

Eq. (18c) from CD98, relative to the temperature, was 
not included in the final network since, being superadia- 
baticity at the bottom of the solar CZ negligible, temper- 
ature itself turned out to be nearly stationary already at 
the beginning of the relaxation. 

As for the diffusive terms Df, the most simple down- 
gradient approximation has been chosen, namely, for each 
generic turbulent quantity Q: 



(3) 



only other non-local treatment ( |Xiong 1985| ) with which 
extensive stellar models have been computed. In fact, in 
Xiong case, an explicit, arbitrary scale length was included, 
and the results almost linearly depended on the value of 
the scale length. After 2 — 3 x 10 6 s, all the six quan- 
tities (including 0) reached a final, stationary distribu- 
tion, clearly showing a thin overshooting region. Figure 
1101 shows the behavior of the turbulent flux J; Figure ITT1 
presents the mean quadratic velocity variance w 2 . Over- 
shooting is present indeed, but its amplitude does not 
overcome 0.02H p , consistently with the solar observational 
constraints. 4 The same relaxations have been performed 
for the CZ of the solar mass track at various evolutionary 
phases, from the first appearence of a radiative nucleus in 
PMS, up to early red giant. The overall features turn out 
to be similar in all cases. It is found that the overshooting 
region is absolutely negligible (< 0.005i/ p ) in PMS during 
the lithium burning phase. This at least ensures us that in- 
clusion of proper overshooting in PMS evolutionary mod- 
els will not worsen the problems with the exceedingly large 
solar lithium depletion shown by todays standard mod- 
els. The thickness of overshooting, then, steadily increases 
during the evolution, reaching ~ 0.06H p at logL/L = 1, 
where the computations have been stopped. In all case, the 
decay of w 2 is very sharp, putting an end to overshooting. 

The only conclusion we can presently draw from these 
first non-local results is that probably the chosen approx- 
imation for the diffusive terms is perhaps not too bad, at 
least as long as the thickness of overshooting is concerned, 
since they are consistent with the observational solar con- 
straints, and do not worsen the lithium problem. 
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Figure 7. Simulations based on the Montalbdn et al. (2003) 
models of different metallicites. A population of 1500 stars 
is randomly chosen between ages of 10 7 and 10 9 j/r defines 
the Bohm Vitense gap at the boundaries indicated by the 
diagonal lines. 




7000 6000 5000 4000 

T 

'et! 

Figure 8. The Pleiades data by Soderblom et al. (1993) 
and Garcia Lopez et al. (1994) (open squares) are com- 
pared with the depletion predicted by the models in Ta- 
ble 2 (full big squares). The models are placed at the 
T e g they would take in an empirical MS, at the Pleiades 
age. Only the upper squares, corresponding to the mod- 
els with a a t m = a- m — 1, are compatible with the data. 
The full line shows the depletion from the models by Ven- 
tura et al. (1998b) computed including the thermal effect 
of a magnetic field on the convective temperature gra- 
dients. The large open square with the error bar repre- 
sents the lithium abundance of the secondary component 
of RXJ 0529.4+0041 (log N(Li)=2.4±0.5, Covmo et al. 
2001). The T e ff at which the point is located (5500K) is 
assumed to be the main sequence T c s of a star of mass 
0.925MV . 
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Figure 9. The figure shows the tracks of 0.8, 0.9, 1.0 and 
I.IMq computed with AH97 atmospheres and aj n = 1.9 
(full lines, left), or a\ u = 1.0 (dash-dotted lines, right). 
The solar location is shown: it is compatible with the solar 
model having a- m — 1.9, but ^jJ_00K hotter than the solar 
model with a- m = 1.0. The location of the secondary com- 
ponent of the binary RXJ 0529.^+00^1 is also shown. Its 
mass is 0.925±0.005Mq , so it is best compatible with the 
a = 1 models 




3x10* 



2x10" 




Figure 11. The mean quadratic velocity variance w 2 . Ve- 
locity decays outside the local Schwarzschild boundary with 
the same slope as inside. When it approaches zero, how- 
ever, the slope becomes steeper, and overshooting comes 
to a full stop. This suggests that all the evolutionary com- 
putations performed up today in the hypotesis of an expo- 
nential decay of the diffusive coefficient outside formally 
convective regions can perhaps (largely?) overestimate the 
amount of extra-mixed matter. 



Figure 10. The non-local turbulent flux J (cgs units) 
around the local Schwarzschild boundary (vertical dashed 
line ) is shown as a function of the logarithm of the pressure 
P (cgs units). The flux is obviously negative in the over- 
shooting region and, contrarily to what one could expect 
according to a thumb rule (perhaps exponential decay), a 
very sharp decay of the flux is found, suddenly ending the 
extra-mixed region less than 0. 02H p below the local bound- 
ary. 



